Effect of spatial bias on the nonequilibrium phase transition in a system of 

coagulating and fragmenting particles 
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We examine the effect of spatial bias on a nonequilibrium system in which masses on a lattice 
evolve through the elementary moves of diffusion, coagulation and fragmentation. When there is no 
preferred directionality in the motion of the masses, the model is known to exhibit a nonequilibrium 
phase transition between two different types of steady states, in all dimensions. We show analytically 
that introducing a preferred direction in the motion of the masses inhibits the occurrence of the 
phase transition in one dimension, in the thermodynamic limit. A finite size system, however, 
continues to show a signature of the original transition, and we characterize the finite size scaling 
implications of this. Our analysis is supported by numerical simulations. In two dimensions, bias is 
shown to be irrelevant. 



PACS numbers: 64.60.-i, 05.70.Ln 



INTRODUCTION 



Systems far from equilibrium can undergo nonequi- 
librium phase transitions between two different types of 
steady states when the parameters of the system are var- 
ied. It is important to know how robust such transitions 
are with respect to changes in the governing dynamics 
and to ask if a signature of the original phases remains, 
even if the transition is lost. The introduction of a spatial 
bias (a preferred direction) is one factor which is known 
to affect the scaling functions and the exponents char- 
acterizing nonequilibrium transitions. Examples where 
bias plays a role include models of extremal dynamics 
jl], 01, the simple exclusion process ||, sand-pile mod- 
els |Q] , directed percolation jl| , interface depinning || , 
reaction-diffusion systems |tJ and random walkers in frac- 
tal media |sj . These examples consist of systems in which 
parameters need to be tuned to reach criticality as well 
as systems that are self-organized critical. In both cases, 
bias either changes the universality class characterizing 
the system ^, g, [| [| ||, @] ° r causes localization H or 
induces boundary driven transitions 

In this paper we study the effect of bias on a re- 
cently introduced model of aggregation and fragmen- 
tation, which was shown j^, |lCJ to exhibit an unusual 
nonequilibrium phase transition belonging to a univer- 
sality class different from models studied earlier || pT| . 
We show that bias in this model plays a different role to 
any of the cases mentioned above in that even an arbi- 
trarily small bias inhibits the phase transition entirely in 
one dimension. Remarkably though, a signature of the 
transition remains and modifies the finite-size behavior of 
the system, and we characterize the scaling implications 
of this. In two and higher dimensions, we show that bias 
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has no effect on the phase transition. 

We define the model on a e?-dimensional hyper-cubic 
lattice with periodic boundary conditions. Starting from 
a random distribution of non- negative integer masses, the 
system evolves in time via the following microscopic pro- 
cesses. In an infinitesimal time dt, (i) with probability 
pdt (qdt), the mass at each site hops to one of its near- 
est neighbors with increasing (decreasing) coordinates, 
(ii) with probability wpdt (wqdt), unit mass is chipped 
off from an already existing mass and added to one of 
its nearest neighbors with increasing (decreasing) coor- 
dinates. Following these moves, the masses at each site 
add up (see Fig. [l]). The dynamics conserves the total 
mass of the system. Hence, the parameters defining the 
system are the bias p — q, the density p and the chipping 
or fragmentation rate w. The case p = q corresponds to 
the zero bias case or the symmetric model while the case 
p 7^ q introduces a preferred direction in the motion of 
the masses and corresponds to the asymmetric model. 

The model may be mapped onto generalizations of 
other well studied models of nonequilibrium statistical 
mechanics. In one dimension the system of masses de- 
scribed above may be equivalently thought of fl(]| as a 
collection of particles and holes on a ring, or a one di- 
mensional interface evolving in time (see Fig. p. In this 
language, the w — > oo limit exactly corresponds to the 
well studied symmetric (asymmetric) exclusion process 
p2| for p = q (p ^= q). In the interface language, this 
corresponds |u| to an interface evolving via the Kardar 
Parisi Zhang or the Edwards- Wilkinson dynamics |l4j for 
the cases p ^ q and p — q respectively. For a finite w, 
the nearest neighbor particle exchanges of the exclusion 
process (or corner flips of the interface growth model) are 
further augmented by long-range moves. 

Models similar to the one studied in this paper have 
also been studied in other contexts. A slightly different 
off-lattice version of this model was studied within the 
rate equation approach in the context of polymer chain 
growth fll5| |. In this case, the mass clusters were thought 
to represent polymers. Although the aggregation of poly- 
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FIG. 1: (a) The model as defined in the text, (b) To obtain 
the hard core lattice model from the masses, lay down the 
masses on their sides, (c) The interface is obtained by replac- 
ing every shaded circle by a line segment in the +45° direction 
and empty circle by line segment in the —45° direction. 
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FIG. 2: P(m) is plotted against m for three values of density. 
The cutoff is seen to depend on p for small values of the den- 
sity while for larger values, it is independent of p. The mass 
in excess of the critical mass forms an aggregate. The solid 
line has an exponent —2. The simulation was done for system 
size V = 500 and w = 3.0. In the inset, the dependence of 
the cut off of the power law on V is shown. The simulations 
were done for p = 15. 



mers and dissociation of single monomers was allowed in 
this model, it lacked the important process of local dif- 
fusion which we include in our study. Various models 
of coagulation and fragmentation, with coagulation rate 
proportional to mass, have been studied in the context 
of gelation |ll| . In these systems, due to the enhanced 
coagulation, a gel (a cluster which contains a finite frac- 
tion of the total mass) forms at a finite time. The ex- 
ponents at the transition point were shown to depend 
on whether the process of fragmentation was present or 
absent. Models of fragmentation in which a fraction of 
existing mass (as opposed to a single particle) may break 



off have also been studied in |17], |18| |19| . In these mod- 
els, it could be deduced that an infinite aggregate that 
contains a finite fraction of the total mass never forms. 
A model with bias very similar in spirit to the model we 
study here, was also studied in the context of traffic flow 
p0| and the phase transition observed numerically was 
interpreted as a traffic jam occurring for large densities 
of cars. 

An understanding of the steady state reached by the 
system may be obtained by considering the limits of only 
diffusion (w = 0) or only chipping (w — oo). Both these 
cases are exactly solvable. In the former case, the system 
maps to the well-studied reaction diffusion system mA + 
nA — > (m + n)A, and the steady state reached is simply 
one in which all the particles accumulate on a single site. 
In the opposite limit (w — > oo), the model is again exactly 
solvable and the system reaches a steady state in which 
the mass is uniformly spread out over the system with 
the probability that a given site has a mass m being 
exponentially distributed. 

The special case p — q, corresponding to zero bias, 
was studied earlier using a mean field approximation M , 



Monte Carlo simulations [|l0| |21{ ], as well as some exact 
analysis |2l]j . It was shown that the system undergoes 
a nonequilibrium phase transition in the p — w plane at 
some critical density p c (w). In particular, P(rn), the 
probability that a randomly chosen site has mass m was 
shown to vary for large m as i) P(m) ~ exp(— m/m*) for 
p < p c , ii) P{m) ~ mT T at p = p c with r = 5/2 and (iii) 
P(m) ~ m~ T + "infinite aggregate" for p > p c , where by 
"infinite aggregate" we mean a single large mass equaling 
a finite fraction of the total mass of the system. That is, 
the mass distribution P(m) changes from an exponential 
distribution to an algebraic one at p c . For p > p Cl the 
mass distribution remains the same as at p c , while all the 
mass in excess of the critical density coagulates together 
forming an infinite aggregate. The mathematical mech- 
anism giving rise to the infinite aggregate was found to 
be very similar to that of equilibrium Bosc-Einstein con- 
densation in an ideal Bose gas. 

Bias is introduced in the motion of the masses on the 
lattice by letting p =/= q. In the limits of only diffusion 
or only chipping, the steady states reached are the same 
as those for the p — q case. Hence, one might expect a 
phase transition in the p — w p lane as before. Further, 
the mean field analysis ||, [l(J does not recognize any 
difference between the two cases and thus predicts the 
same behavior as in the zero-bias case. Earlier Monte 
Carlo simulations of this system in one dimension did 
indeed seem to suggest the existence of a phase transition 
similar to the symmetric case, though with r « 2.0 Ml. 
The results of these simulations are summarized in FigT|2|. 
For a fixed lattice size and w, when the total mass in the 
system is increased beyond a certain critical mass, the 
formation of an aggregate is observed. The power law 



3 



regime has a lattice size dependent cutoff (see inset of 
Fig. ||) as in the symmetric case. However, the value of 
r (~ 2.0) is a puzzle since a finite system density implies 
that r should be strictly greater than 2 (the first moment, 
which is the density of the system, would diverge if r < 
2). 

In this paper, we analyze the model for p =/= q and show 
that the apparent existence of a phase transition in the 
Monte Carlo simulations is purely a finite-size effect in 
one dimension, with the critical density p c (V) for fixed 
w diverging with system size V as log(F). We argue that 
the exponent r is exactly 2. In two and higher dimensions 
however, a transition from an exponential to an aggregate 
phase does occur at finite critical density. 

The rest of the paper is organized as follows. In Sec. ||, 
we show analytically that an aggregate phase cannot ex- 
ist in the one dimensional system, thereby showing that 



III 



there is no phase transition in one dimension. Sec. 
contains numerical evidence for the results of Sec 
from Monte Carlo simulations of the model. In Sec 
the exponents describing the probability distribution for 
mass in one dimension are characterized. Section [v| con- 
tains a numerical study of the two dimensional problem. 
Section VI contains a summary and conclusions. 



II. ARGUMENTS FOR NO PHASE 
TRANSITION IN THE PRESENCE OF SPATIAL 
BIAS IN ONE DIMENSION 



In this section, we prove that in the presence of a spa- 
tial bias, an aggregate phase cannot be present at any 
finite density in one dimension. We do so by assuming 
that an aggregate phase exists, and then showing that 
this leads to certain contradictions. To proceed, assume 
that a single infinite aggregate exists in the system and 
that the rest of the system is at a finite critical density 
p c (analogous to the symmetric case). Consider now a 
frame of reference that is attached to this aggregate. Let 
r) k (t) be the mass transferred in an infinitesimal time dt 
at time t from a site k lattice sites away from the aggre- 
gate to a site k ± 1 sites from the aggregate. Then 

m k (t+dt) = m k (t)+r ] +_ 1 +i ll : +1 -r]+-r]' k +a++a^, (1) 

where m k (t) is the mass at a site k lattice sites from the 
aggregate at time t, and a k is the change in mass due 
to the diffusion of the aggregate. The time dependence 
of the variables on the right hand side of Eq. ([!]) have 
been suppressed for the sake of clarity. Equation ([!]) 
accounts for all the ways in which the mass on a site 
k sites away from the aggregate can change under the 
dynamics. These changes are either caused by the mass 
transfer to and from neighboring sites as exemplified by 
the rjk's or by the motion of the aggregate, which leads 
to a relabeling of sites, as depicted by the a k s. 



From the definition of the model, it is clear that 

' m k with prob pdt, 

with prob 
otherwise, 



Vk = { 1 - <W,0 with prob pwdt, (2) 



and similarly for rfj^ , with p replaced by q in Eq. (||). 
When the aggregate hops, the sites have to be relabeled 
and this leads to 



m k+1 - m k with prob pdt, ,„.. 
otherwise, ^ ' 



and correspondingly for , with k + 1 replaced by k — 1 
and p by q in Eq. (Bh . Taking averages on both sides in 
Eq. (|l|) and setting the time derivatives to zero in the 
steady state, we obtain 



dpk 
dt 



= w [psfc_i - (p + q)s k + qs k+ i] + {p + q)x 



[pk-i - 2pk + Pk+i + Po(2fa,Q - fa-i - fa,i)] ■ (4) 

where s = 1, and p k = (m k ) and s k = (I - S mkfi ) 
are the average density and occupation probability of a 
site k lattice sites away from the aggregate respectively. 
A point to note about Eq. (^) is that bias plays a role 
only in the terms coupling to the chipping rate w. It 
can be checked that introducing a bias only in the dif- 
fusion move and keeping the chipping move symmetric 
does not change the behavior of the model from the fully 
symmetric version p = q. 

We now consider Eq. (|j) in the steady state when the 
time derivative is set to zero. The set of linear equations 
in Eq. (Q) may be solved for on a finite or an infinite 
lattice to obtain the p k s in terms of the s k s. In the 
former closed form expression for the densities is 

easily obtained. However, it is more informative to look 
at the equations for an infinite lattice, where the sites 
with negative (left of aggregate) and positive (right of 
aggregate) indices may be treated separately. For this 
case, we obtain for n > 0, 



+ ^ — :( ( ? s i -p)^ : — s - 



p+q- p+q 
W (P ~ l) „r\ i W P 



p + q 

P-i - 

w{p - q) 
p + q 



X>k-*c) + --X-(l-*»), (5) 

^ p + q 



k=l 
W 



p + q 



w{p-q) l 
-iPS-i -q) : si 



p + q 



fc=-i 



wq 
p + q 



(l-s-»),(6) 



where s r r = linin^oo s n and s 



lim r . 



s — j 



By assumption, we require that p n tends to a finite 
value as n — > oo. For this, we require first that the 
term proportional to n in Eqs. (^]) and ([|) vanishes and 
secondly that s± n approaches its asymptotic value faster 
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than 1/n. The first condition expresses p\ and p-\ in 
terms of S\ and s_i as 



Pi 



P-i = 



w(j> - qsi) - w(p - q)s r c 
p + q 

w(q — ps-i) — w(q — p)t 
p + q 



(7) 
(8) 



which when substituted in the k — equation of Eq. (^) 
leads to 



(9) 



This is consistent with the fact that far from the ag- 
gregate, the occupation probability is the same on ei- 
ther side. The second condition, that s± n approaches its 
asymptotic value faster than 1/n, will be useful for de- 
termining the exponents and will be studied numerically 
in Sec. [ffl. 

A important point to note from Eq. (||) is that the two 
cases p = q and p ^ q are quite different. For p = q, the 



cumulative sum over the s„'s vanishes and p n depends 
only on the site occupancy of the site n. However when 
p ^ q, the sum plays a role in determining the value of 
the site density and it becomes important to understand 
the behavior of s± n as a function of n. We do that in 
Sec. |m. 



We now examine the two point correlations in the pres- 
ence of the aggregate in order to obtain further relations 
between the pk's and Sfc's. Using these relations we will 
be able to show a contradiction. Consider the two point 
correlations in the aggregate frame of reference. Let 



C r ,k = (m r m r+k ). 



(10) 



The equations governing the temporal evolution of the 
two point mass-mass correlations may be derived by con- 
sidering the mass transfer between two neighboring sites 
(see Eq. (Q)). Multiplying together Eq. (|l|) for m r and 
the corresponding one for m r+ k, keeping terms up to or- 
der dt, and taking averages, we obtain 



dC r ,k 
dt 



dt 



pC r+ i tk + qC r+li k-l + pC r ,k-i + qCr.k+l + pCV-i.fe+l + qC r -i ik - 3(p + q)C r ^ k 
+ w [-pD r -x <k+ i + {p + q)D rt k - qD r+lik -i - pE r ^-i + (p + q)E r>k - qE r ^+i] 

- 4,1 [qCr+l.O + pCrfi + Wps r + WqS r+ i] , T, k = 1, 2, . . . , 

(p + q) [CV+i,o - 2CV,o + C r _i, ] + 2pCV-i,i + 2gC r ,i - w [2pD r _ 1A + 2qE r>1 ] 
+ w[(p + q)s r + ps r _i + qs r+1 ] r = 1, 2, . . . , 

I 



(11) 
(12) 



where 



D r ,k = (S mr ,om r+k ) fc=l,2, . 
E r ,k = {m r d mr+k fi) ft = 1,2,. 



(13) 
(14) 



Also, for Eqs. (11) and (12) to be valid when r — 1, we 
need to set Co,fe = A),fc = 0. In the steady state, the time 
derivatives can be set to zero. Summing Eq. ( |TT| ) over all 
r, k = 1,2,..., and subtracting from this 1/2 times the 
sum over r = 1, 2, ... in Eq. (|l2|), we obtain 

oo oo 

(p + q)C lfi = wp-2(p + q)^Ci t k + 2wq^D ltk + wqs 1 . 

k=l k=l 

(15) 

The left hand side of Eq. ( ]l5| ) is a finite number while 
the right hand side has two infinite sums. These two 
sums can add up to give a finite value only if asymptot- 
ically the terms in the two summation are equal. Using 
the fact that the two point correlations decouple when 
the separation between the two points become large, we 
obtain 



Equation ( pf ) and (0) have to be simultaneously satisfied. 
This is possible only if 



w(p -q)(l- s c ) = 0. 



(17) 



(p + q) pi = w <?(i ~ si). 



(16) 



The fragmenting rate w being equal to zero is the trivial 
limit in which the entire mass in the system coagulates 
together to form an infinite aggregate. When there is a 
bias p ^ q, then for finite w, the only way Eq. ( p~7[ ) can 
be satisfied is if s c = 1 ; i.e, the occupation probability 
far away from the aggregate is 1. This can occur only 
if p c = oo, which contradicts our initial assumption that 
p c is finite. Another way of seeing a contradiction is 
to consider Eq. (||) when n — > oo. Setting s c — 1, we 
obtain that the densities far away from the aggregate 
become negative for p > q. Thus, for any finite w and 
p we always have a contradiction and hence our initial 
assumption of an aggregate existing at finite density is 
proved wrong. This proves that for any finite p and w, 
an aggregate phase does not exist in one dimension in 
the thermodynamic limit. For the symmetric problem 
(p = q), there is no contradiction between Eq. (16) and 
(^), since Eq. ( |l7j ) is automatically satisfied. 
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FIG. 3: /o c (V) diverges with V as \og(V). The simulations 
were done for the fully asymmetric model p = 1 and q — 0. 



III. NUMERICAL CHECKS IN ONE 
DIMENSION 



The results of the previous section thus show that an 
infinite aggregate cannot exist in an infinite system when 
there is a non zero bias. However, numerical simulations 
of finite size systems do point to the existence of an ag- 
gregate phase (see Fig. ||). Nevertheless, there is no con- 
tradiction with the results of Sec. [Oj, provided the critical 
density p c (V) diverges with the system size V. In this 
section, we study numerically the system size dependence 
of p c (V) in one dimension. 

In earlier studies jlO|, the system size dependence of 
p c (V) was not investigated because there was no system- 
atic way of making an accurate numerical measurement 
of the critical density. To measure p c (V), we adopted the 
following procedure. For a fixed lattice size, we start the 
system with a density much higher than that required to 
form an aggregate. The system is then allowed to reach 
the steady state and the biggest cluster is identified as 
the aggregate. We then measure the density in the rest 
of the system (excluding the aggregate) and use the fact 
that the state of the rest of the system resembles that 
at criticality. In Fig. g the system size dependence of 
p c (V) is shown on a semi-log scale. From the numerical 
evidence, we conclude that p c (V) ~ log(V). 

As a further check, we study the occupation prob- 
ability numerically. From Eq. (|B|), the dependence of 
p n on XX s ±fc — s c)i taken together with the fact that 
Pc{V) ~ log(V) implies that 



\S±k 



a± 
|fc|*' 



k > 1, 



(18) 



with x = 1 which in turn implies than s c (V) converges 
to its asymptotic value as 1/V. Both of these require- 
ments are consistent with numerical simulations (see 
Fig. ^). The simulations were done for the fully asym- 
metric model p = 1 and q = 0. 
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FIG. 4: Sfc converges to s c (~ 0.2775 in this case) as 1/k. The 
simulations were done for lattice size 4000 and w = 1.0. The 
inset shows the finite size correction to s c . 




FIG. 5: The probability distribution for mass for different 
values of p. The inset shows the scaling plots when the dis- 
tributions are scaled as in Eq. (po|). 



IV. 



PROBABILITY DISTRIBUTION 



From the analytical and numerical evidence of sec- 
tions H and III it is clear that the system is sensitive to 
the manner in which the limits M — > oo and V — > oo are 
taken, where M is the total mass in the system. When M 
is increased beyond M c (= p c V) keeping V fixed, an in- 
finite aggregate does form in the system. In this regime, 
in analogy with the symmetric problem |2l]] , we can then 
write a scaling form for the probability distribution for 
M > M c as 

P(rn, V) = ±f (£L) + I*(m - (M - M c )). (19) 

Since the mean mass in the power law part of the dis- 
tribution scales as \og(V) (see Fig. ||), we immediately 
derive r — 2. In addition, from the consideration that 
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there is only one aggregate, the two exponents r and <f> 
are known to obey the scaling relation 0(r — 1) = 1 
This implies that = 1. The fact that the cutoff of the 
power-law distribution scales as V, and not as a smaller 
power of V as in the symmetric case 



21 



is consistent 

with the fact that the transition does not exist for large 
V. 

What is the behavior of the system when the order 
of limits is reversed and M, V — > oo keeping p = M/V 
fixed? In this case, we make the reasonable ansatz, 



1 / to 

hm P(m, p) = —g — 



(20) 



0.5664-p c (V) 
0.2865/x**0.25 



100 



1000 



with the same exponent r = 2. The requirement that 
(to) = p taken together with Eq. (|2^) implies that the 
cutoff to* ~ e ap . Scaling plots of the probability distri- 
bution for various values of p scaled as in Eq. (^0|) are 
shown in Fig. |5|. 

From simulations, it is seen that the function g(x) ~ 
const when x — > (See inset of Fig. ||). This taken to- 
gether with the exponential divergence of the mass cut-off 
implies that for infinite p, P(m) is a pure power-law. This 
is similar to the steady state of the Takayasu model p^] 
for river networks where mass aggregates in the presence 
of a constant influx of particles, and the steady state has 
a nontrivial power law distribution. 



NUMERICAL SIMULATIONS IN TWO 
DIMENSIONS 



FIG. 6: The convergence of the critical density to a finite 
value with increasing system size in two dimensions is shown. 
The simulations were done for the fully asymmetric model 
p = 1, q = and the fragmenting rate w = 1.0. 
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The arguments used in Sec. |lj to prove that there is no 
phase transition in one dimension are very specific to one 
dimension and cannot be extended to two and higher di- 
mensions. Instead, in this section, we numerically study 
the model in two dimensions. First, we measure the crit- 
ical density p c (V) for lattice sizes varying from 8x8 
to 32 x 32 using the same method as that used for the 
one dimensional simulations. In this case it is seen that 
p c iy) does converge to a finite value when the system size 
is extrapolated to infinity (see Fig. ^). Hence, a phase 
transition does exist in the infinite system limit. 

We now show that for the system with spatial bias 
the value of the power law exponent t in two dimensions 
is close to the zero bias case. It is difficult to make a 
direct measurement of r because the cutoff to the power 
law grows slowly with the system size and is not large 
enough ((cutoff ~ 200 when the system size ~ 2000) for 
an accurate measurement of the slope. For instance in 
the zero bias case direct measurements of r || [To) gave a 
value close to 2.3 while indirect numerical methods 
showed that r is close to 5/2 in all dimensions. Therefore, 
rather than measure r directly for the asymmetric model, 
we compare the simulation results (see Fig. 0) for the 
fully asymmetric problem with those for the symmetric 
problem with the same parameters. The slopes of the 
two curves are comparable. Hence, we conclude that the 
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FIG. 7: The value of the power law exponent r for the fully 
asymmetric model in two dimensions is very close to that 
of the symmetric model in two dimensions. The simulations 
were done for a 30 x 30 lattice with w = 1.0 and p — 10.0. 

exponent r for the asymmetric model in two dimensions 
is very close to 5/2, as in the symmetric model. 



VI. SUMMARY AND CONCLUSIONS 

In summary, we have investigated in detail the effect 
of introducing a spatial bias on the nonequilibrium phase 
transition in a model of coagulation and fragmentation. 
We show analytically that the phase transition is inhib- 
ited in one dimension. However a signature of the two 
original phases remains and the scaling implications of 
this are characterized. We have also resolved the puzzle 
of the exponent r being very close to 2. In two dimen- 
sions, the phase transition is shown numerically to exist. 

We now give a more intuitive explanation of why the 
phase transition gets curbed in one dimension but not in 
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two dimensions. In this model, there are two competing 
processes. While the diffusion move creates larger and 
larger masses by coagulation, the fragmentation move 
tends to create smaller masses as well as inhibit the for- 
mation of large masses. If the diffusion move was to be 
considered by itself, then a cluster of size I would be cre- 
ated in time of order I 2 . In one dimension, if the fragmen- 
tation move was to be considered by itself, then a fluctu- 
ation of order I would be dissipated in time of the order 
? 3 / 2 for the asymmetric model and of the order I 2 for the 
symmetric model. These exponents are known exactly 
because of the existing exact analogy in one dimension, 
between only fragmentation with (without) bias and the 
asymmetric (symmetric) exclusion process. Thus, for the 
asymmetric problem, fragmentation always wins out over 
diffusion and we only have an exponential phase. How- 
ever, in two dimensions, bias is irrelevant for the frag- 
mentation move and hence a fluctuation of order I gets 
dissipated in time of order I 2 , which is of the same or- 
der as the time required to create a cluster of size I by 
diffusion. 

To carry this argument further, we can study the sym- 
metric problem by slowing down the diffusion rate. This 
can be done by introducing a mass dependent diffusion 
rate ~ m~ a with a > 0. The above arguments would 



then imply that this dynamics ought not to have a phase 
transition for any a > 0. This is indeed the case and it 
can be shown that the phase transition does get curbed 
in all dimensions Q, as predicted. 

There remain several interesting questions to investi- 
gate further. While, we have numerically shown that 
t = 2.0, it would be interesting if it could be derived from 
first principles by solving the model. Further, for an in- 
finite system, the probability distribution of the masses 
has the form (Eq. @) 

P{m) ~ -i- e -^ me "° P , m»l, 

TTl 

where (3 could depend on p and w while a depends only 
on w. The origin of the length scale e ap under this dy- 
namics is an interesting point that remains to be under- 
stood. Also, since the phase transition in one dimension 
is a finite size effect, the implications of the traffic jam 
that was seen in |^0j need to be reexamined. 
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